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Abstract. In this paper we present a stability criterion for finite measure- 
valued stochastic recursions, generalizing Loynes's Theorem to spaces of mea- 
sures. This result provides conditions for the reach of a "total stationary 
state" for the queue with an infinity of servers and the single-server SRPT 
queue. Indeed, we give in both cases a condition of existence of a stationary 
measure- valued recursive sequence characterizing the queueing system exhaus- 
tively. 



1. Introduction 

In this paper, we address the question of stationarity in the general ergodic frame- 
work for queues, that admits a representation by a measure- valued stochastically 
recursive sequence (SRS). We axamine the good topological conditions for Loynes' 
backward recurrence scheme to handle the question of existence of a stationary ver- 
sion (in the sense Loynes [2 J puts it) for such an SRS. After some preliminaries, we 
prove that these conditions (essentially, that all increasing sequence in some sense 
converges) are met when the space of finite non-negative measures is endowed with 
the weak topology and the "strong" partial ordering of Stoyan. 

Then, Loynes' construction allows one to prove the "total" stability of a given 
system whenever a measure valued SRS describing it exhaustively enjoys some 
kind of monotonicity. Hereafter we follow this construction for two queueing sys- 
tems : the queue with an infinity of servers and the Shortest Remaining Processing 
Time (SRPT) single server queue. Some applications of these results, such as weak 
convergence of the congestion and workload sequences of these queues (which are 
projections of the measure- valued SRS for adequate test functions), are given. 

This paper is organized as follows. After some preliminaries in section we 
state our Loynes' type theorem for measure- valued SRS in section 03 In section [4] 
we study a particular real-valued stochastic recursion, that will be usefull in the 
sequel. In section 03 (respectively, section 03), we study the stability of the queue 
with an infinity of servers (resp., of the SRPT queue) through its measure- valued 
representation. 



2. Preliminaries 

Let Mj and Cb denote respectively the set of positive finite measures on R+ 
and the set of bounded continuous functions from R to R. Equipped with the 
weak topology tr ( Mt,Cb ) , Mt is Polish (see [1]). Let £ be the zero measure on 
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K (i.e., such that C(^) = for any Borel set «8 on R). For any /i G M+ and 
any measurable / : M — > R, we classically write (fi, /) := f f dfi. We denote 
C (Mj^ , the space of continuous mappings from My into itself. Let the set My 
be endowed with the increasing partial integral order ^ : for any two /x, v G My , 
M ^ v if (/x, /) < ( v i f) f° r an Y measurable non-decreasing function / such that 
these integrals exist. Remark, that ( ^ /x for any /x G My . Let us denote for any 
y G R and any measurable / : R — > R, r y /(.) = /(. — > !/ }. Then, for any 
it G My , r^zx denotes the only element of My s.t. (r^it, /) = (fi,T y f). 

Let now Al G My be the subset of finite counting measures on R+. Any 

/x G A4\{(} reads /x = X^f=i + ^ ^a.O)> where is the Dirac measure at a; G R+ and 

aiO) < a 2 (/i) < ... < a^irpOx). Then, r y (/i) = ^(aO-i/^O*)^} and for 

any two /x, ^ G .M\ {£}, fi ^ f whenever 

(i) mr;) < f(R+), 

(m) for all i = 0, ...,/z(R+) - 1, a M ( R ^)_i(^) < a^(M;)-i(^)- 

Let us finally denote for any /x G -M\{C}, to(zx) = ai(/x) and M(/i) = a M (Rp(/x), 
the smallest and the largest atom of /i, respectively. 

Lemma 1. j4njy sequence o/M| iftai is ^-increasing and bounded above converges. 

Proof. Let {xt,i}„ gN be a ^-increasing sequence of My that is bounded above 
by it G My . Then, as easily seen the sequence of non-increasing real functions 
{[j, n ([., oo))} n£N tends pointwise, and hence (this is Diniz Theorem), uniformly to 
a non-increasing real function / that is right continuous and has a coutable number 
of discontinuities. Moreover /(0) < /x(R+) < oo, and we can fully characterize 
a measure xt* G My setting xt* ((0,x)) = /(0) — f(x) for all x G W + . In particu- 
lar, sup^gR. |/x™ ((0, x)) — ii* ((0, a;)) | — ► 0, hence jjb n tends to /x* in total variation, 
which completes the proof. □ 

3. Stationarity of measure- valued stochastic recursions 

Consider a probability space (O, T, P°) , embedded with the measurable bijective 
flow 9 (denote 6* _1 , its measurable inverse). Suppose that P° is stationary and 
ergodic under 9, i.e. for all 51 £ J, P° [tf^Sl] = P° [21] and 021 = 21 implies 
P° [21] = or 1. Note that according to these axioms, all ^-contracting event 
(such that P° [2t c n^-^] = 0) is 0-invariant. We denote for all n G N, 9 n = 
9 o 9 o ... o 9 and 9~ n = 9~ 1 o 9~ l o ... o 9~ 1 , and say that two sequences of r.v. 
couple when there exists a P°-a.s. finite rank N such that they coincide for all 
n > N. 

Let $ be a C(My )-valued r.v., k be an My-valued random variable (r.v. for 
short), and define the following sequence of Mt-valued r.v.. 



r Mo = 

\ it£ = * o 0" (z^) , „ G N. 
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The existence of a stationary (i.e., compatible with 6) version of {^n} n <£N amounts 
to that of a solution of the following equation in the space of Mjr-valued r.v.. 

(1) ko0 = $(k). 

We have the following result. 

Theorem 1. The equation (Qp admits a solution whenever $ is a.s. <-non de- 
creasing and the sequence °^~ n }„ e N * s a - s - di-bounded. 

Proof. With Lemma Q] in hands, we can construct a solution to {1} following the 
classical Loynes's backwards scheme ( [2], [4] p. 107): the sequence o &~ n } neN 
is P°-a.s. non decreasing and bounded, and thus converges a.s. to the ^-minimal 
solution of {J), which belongs P°-a.s. to Ml". □ 

4. Study of a real valued stochastic recursion 

Under the foregoing definitions and assumptions, let a > 0, (5 > be two 
R+- valued integrable random variables (r.v. for short), and the sequences of r.v. 



Let Z be an a.s. finite 



-valued r.v.. Con- 



(E) 



sider the following stochastically recursive sequence : 
Y * = Z, 

Y%+\ = [max {r r f , a n } - j3 n ] + for all n £ N. 

In this section we are interested in the existence and uniqueness of a stationary 
solution } n£N to the recursion (E), in that we are looking for an a.s. finite r.v. 
Y such that for all n G N 

r„ r = yof. 

In particular, Y would be a solution to the equation 
(2) Yo6= [max {Y, a}- f3] + . 

We have the following: 

Theorem 2. There exists an only P°-a.s. finite solution Y of {!)), given by 



(3) 



Y 



sup 



J 

E 

i=l 



P-i 



Moreover, for all P°-a.s. finite and nonnegative r.v. Z , the sequence {^^} neN 
couples with {Y o ^™}„ gN . 

Proof. Existence and finiteness 

Equation |f2]) can be handled by Loynes's construction ( [2]) since the functional 
[max{., a} — f3} + o6* _1 is P°-a.s. continuous and nondecreasing. Loynes's sequence, 
defined by {^°6 ,_ "} neN , 

is given for all n > 1 by 



max ct-j 

3 = 1 \ 



for all n G N. 



Thus, a solution to (fSJ is given by 

Moo := lim M n = 



sup 

jGN* 



3 
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and for any other solution Z, 

(4) P°[Z>Moo] = l. 

Furthermore, Birkhoff's ergodic theorem clearly entails that Moo is P°-a.s. finite. 

Uniqueness 
We have 

(5) P° [Z < a] > 0. 

Indeed, if P° [Z = 0] > 0, there is nothing to prove. Now if P° [Z = 0] = 0, 
P° [Z < a] = entails that on a P°-almost sure event, 

Zo9>0o [max {Z,a}-/3] + >0«Zofl = max {Z,a} - (3 = Z - (3. 

We thus have 

E° [Z o 9 - Z] = -E° [/3] < 0, 

which is absurd in view of the ergodic lemma ( [4], Lemma 2.2.1.). Thus |(5]) is 
verified. Since 

{Z < a} = iz o 9 = [a - 0\ + \ C {Z o 9 < o 9} , 
JH) implies that 

(6) P° [Z < Moo] = P° ° < Moo ° 0} > 0. 
But on {Z < Moo}, 

Z o 9 = [max {Z, a} - (3} + < [max {Moo, a} - (3} + = o 9, 
which means that the event {Z < M M } is ^-contracting. Thus, with iJH), 

P° [Z< Moo] = 1, 

and with (gj), 

P° [Z = Moo] = 1. 

Coupling 

The developped form of Y® is given for all n G N by 



(7) ^ = 



n-l 

max | a 

3=0 



n-l 

3 



«=3 



with the conventions max^ . = and X)j=o • = 0- < Z < oo, P°-a.s..We 
have for all n £ N 



(8) Y* = max \Z- F ? 



n f 7 



i=0 



as easily seen by induction. Indeed, 

Y Z = Z = max{Z,0}, 
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and if (El) holds for some n£N, 



Y z 

1 n+l 



n-1 

max { Z — j3- 

i=0 




Thus, 

{y n z ^ y° for all n G N} = {F,f > F° for all n G N} 



Y r f = Z -J2( 3 i>° f or all n G N 



k i=0 

But the latter event is of probability zero since it contradicts 

n-l \ ~> 



lim - U-Vft 

\ i=0 



E° < [> , 



which is P°-almost sure from Birkhoff 's theorem. Hence, since two sequences driven 
by the same recursive equation couple from the instant when they coincide for the 
first time, there exists a random indice 

t z -° = inf {k G N, Y r z = Y° for all n > k} < oo, P°-a.s.. 

This is true for any finite and nonnegative r.v. Z, and in particular for Z = 



Y. Thus, for all rank larger than r 



Z,Y 



{t*.V.°}<oo, {Y n z \ 



and 
□ 



{Y r Y = Y o n } neN coincide: these two sequences couple. 

Corollary 1. For all finite and nonnegative r.v. Z, there exists P°-a.s. an infinity 
of indices such that Y z = if and only if : 

sup I a-j - ^2 P-i I - > °> 



(9) 

_ \ 2=1 
Proof. In view of the previous coupling property, 



1 " 

i=l 

1 ( n t z - y -1 t z y -1 

= „ 1 ™ ~ E 1 { y ° es =°}+ zZ 1 {y i z =o}- zZ 1 {r°(>*=o 



P° [Y = 0] , P°-a.s. 



Hence 



E 1 {v i *=o} = TC < P °- a - s - ^= P ° t y = °1 > °> 



which is condition J9]). Now if ([9]) does not hold, on a P°-almost sure event, 
the sequence j^rf } ngN coincides after a certain finite rank with the R+*-valued 
sequence {Y o #™} ngN , and hence never reaches after this instant. This event is 
9ti -invariant, and thus P-almost sure. □ 
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5. Functional stability of the G/G/oo queue 

Let P,0t) be a probability space furnished with a bijective flow (9 t ) t>0 , 

under which P is stationary and ergodic. Define on fl the ^-compatible simple 
point process (N t ) t>0 of points ... < T_2 < TLj < To < < T\ < T2 < marked 
by a sequence {c n } neZ . Also denote for all n £ Z, („ = T n+ i — T n , and suppose 
that the generic r.v. a and £ are integrable. Consider the G/G/oo queue fed by 
N(a): there is an infinity of servers, thus e.g., the customer C n arrived at time T n 
is immediately attended and spend in the system a time equal to the duration a n of 
his service. Assume that C n is present in the system at time t, i.e., T n < t < T n +a n . 
The remaining processing time (i.e., time before the service completion) of C n at t 
is hence given by T n + a n — t. The profile of the system at t is the collection of the 
remaining processing times of all the customers in the system at t, and can thus be 
represented by the following measure, which P°-a.s. belongs to M. 

£ if the system is empty at t, 

x{t) 

53 <W) ifnot > 

i=l 

where < m = oc\{t) < a^t) < ... < ctx(t)(t) — M (n(t)) are the remaining 

service times of the X(t) customers present in the system at t. Then (l*i(t), t > 0) 
provides an exhaustive representation of the system in that it keeps track of all 
present customers at a given time. Note, that a diffusion approximation of a scaled 
sequence of such profile processes under Poissonian assumptions is presented in [5]. 
In particular, for any t > 0, the congestion X(t) (i.e., the number of customers 
in the system at time t) and the workload W{t) (i.e., the quantity of work in the 
system at time t, in time units) read 

f X(t) =</«(*), 1), 

I W(t) =Eii' ) a»(*)^(t) = <M*U>, 

where / is the identity function. The processes (fi(t),t > 0), (X(t),t>0) and 
(W(t), t > 0) have rcll paths and denoting for any t > 0, /i(i— ) = lim s f| t /i(s) (and 
accordingly, X(t— ) and W(t— )), we have for all n G Z 

(10) ^(*) = Tt-T n 0* (Tn-) +S a J, te [T n , T n+1 ) , P - a.S. 

Let (CI, P°, 0) be the Palm space of N, Denote 9 := 9t x and for all n G Z, 
0™ := # o 6* o ... o 9. Denoting £ := £0 and cr := ooj we have for all n G Z, := £ o 0™ 
and cr„ := cr o n . For all Mj^-valued r.v. k, denote for all n G N, n n = (i (T n — ) 
provided that /z K (To— ) = k (and accordingly, A^ and for any integer- valued 
r.v. A and any real- valued r.v. W). 

The problem of finding a stationary version of the process (//(t)) t>0 amounts 
to that of finding a sequence {/x£} ngN compatible with the discrete flow 9, which 
amounts in turn in view of (flu]) to finding a P°-a.s. finite random variable n of Mj^ 
such that 

(11) Ko9 = Tt.(n + 5 a ). 
We have the 



tigZ 



i+ff„-i l{T„<t<T„+o-„} 
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Theorem 3. The equation f2P admits a finite solution, given by 

oo 

Moo = 5 (D-i-T,Ui t-i) 1 {n-i>Y,U t-j}' 
t=i 

Moreover, provided that 
(12) P' 
holds, this solution is unique. 

Proof. It is a straightforward consequence of Birkhoff's ergodic theorem that 
P° [Moo eM] = P' 

This ^-contracting event is thus P°-almost sure. The first statement then follows 
from Theorem [TJ the mapping fi \—> (fi + 8 a ) is P°-a.s. continuous and non- 
decreasing, and for all n G N, 

n 

Thus /ioo is the -^-smallest upper-bound of {/in} ngN and hence the -^-minimal 
solution of ((TT)l . Let now re be a solution of ifTTj) . It is easily seen that aAr( K )(re) 
satisfies the stationary equation 

«jv(«)(k) °# = [max{a w(K) (re),cr} + , 

which is of the type ((2|)- Hence (Theorem El uniqueness result), ajv(K)( K ) = 
a N{^ aa ) (Moo) • Since /ioo is a -<-minimal solution, 

{re = /ioo} _> {« = 0} = {aAr( K )(re) = 0} = {a N ^) (/ioo) = 0} , 

which has a strictly positive probability whenever (fl2|l holds. The uniqueness is 
proven in this case since {re = /ioo} is ^-contracting. □ 

Corollary 2. TTie following convergences in distribution hold: 
(i) M° — > Moo, 

(a; x° Moo (k;) , 

(»«,) W° — » J xd/j,oc(x). 

Proof. The convergence (i) classicaly follows from Loynes's construction, since for 
all bounded continuous function F : M. — > M: 

E° [F ( M °)] = E° [F (m° o 9-)} — E° [F ( Moo )] . 

Now, the mapping /i i— > m0^+) is continuous from M to M provided that /i has no 
atom at 0, which is the case. This continuity argument holds for (iii) as well since 
the sequence (Wi!} n£N is tight. We therefore can apply the Continuous Mapping 
Theorem (see [1]), to deduce (ii) and (iii) from (i). □ 



sup 



j'=i 



> 



G PT.D-i-j^C-j- > I < oo 



> 0. 
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6. The SRPT queue 

Under the settings and notations of the previous section, consider now a G/G/l 
queue fed by the input N(a). We assume furthermore that the single server obeys 
the Shortest Remaining Processing Time (SRPT) service discipline, i.e., it always 
gives a preemptive priority to the customer having the least remaining service 
time. Let (v(t),t > 0) denote the profile process of the queue, keeping track of 
the remaining service times of the customers in the system at current time. This 
process is governed by the following dynamics: for all n S Z, 

(13) v(t) = T t _ T „ (v{T n -) + <5 CT J , te[T n ,T n+1 ), P - a.s. 
where for any x e M.+ and any fx e M, 

t,(m)= g r {(,-Er.W(A0) + }^' 

setting J2°k=i ■ = 0- Hence there exists a stationary version of {^n}„ eN provided 
that for some .M-valued r.v. n, 

(14) K0d=T^(K + S a ). 

Theorem 4. If E° [cr] < E° [£], (f!^ arfmiis a solution v x S X, P°-a.s. 

Proof. Denote as above {^n}„ S N, the profile sequence provided that Uq = k. The 
result is a consequence of Theorem Q] since the mapping /i i— > (/i + 5^) is P°-a.s. 
^-non decreasing and continuous, and since {y^ o #~"} neN is ^-unbounded only if 
the congestion sequence or the workload sequence initiated at is <-unbounded, 
which is of probability in view of Loynes' theorem for G/G/l queues. □ 
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